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Abstract 

We work in a chiral invariant quark model, with a condensed vacuum, char- 
acterized by only one parameter. Bound state equations for the nucleon and 
A are solved in order to obtain an updated value of their radii and masses. 
Nucleon-nucleon S-Wave scattering is studied in the RGM framework both 
for isospin T=l and T=0. The phase shifts are calculated and an equivalent 
local potential, which is consistent with K-N scattering, is derived. The re- 
sult is a reasonable microscopic short range repulsion in the nucleon-nucleon 
interaction. 
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I. INTRODUCTION 



The NN central repulsion is a prominent experimental feature of NN scattering and 
through the years it has been the subject of many theoretical papers. The explanation to 
why this NN repulsion happens falls into two main classes: the skyrmion picture and 
the resonating group method approach (RGM) |2|-|J. Both approaches yield a sizeable NN 
central repulsion and the question arises whether or not these models can also explain other 
aspects of the NN force and whether and if they can be brought to describe other hadronic 
phenomena. This extensive program is not completely met yet up to the present day . 

An effective model |]5]|| has been developped in a series of papers, with just one scale 
besides the quark mass parameter, which amounts to a BCS theory of "Cooper-like" quark 
antiquark pairs with the quantum numbers of the vacuum: 3 P . This model embodies in 
a natural way the physics of spontaneous chiral symmetry breaking SxSB and therefore 
has a "correct" chiral limit so that it is expected to fare well in describing low energy 
phenomenology. In our model the mechanism of SxSB is triggered by the effective quark 
potential whose strength besides governing the extent of the quark condensation also fur- 
nishes a natural scale (the typical hadronic size) which turns out to be crucial to understand 
exotic scattering cross sections [^]. In this way we were able to describe a considerable 
number of experimental results ranging from hadronic spectrum and exotic scattering cross 
sections to the coupling of pions to nucleons and deltas ||. As for the main weak point 
of this model we must refer that it yields a too small f n . This is a consequence of lack of 
covariance [JToJ] together with the absence of coupled channels in the BCS treatment of this 
many-body problem. For recent developments on spontaneous chiral symmetry breaking in 
quark models, see ||. 

So it is interesting to see whether we can obtain a reasonable NN repulsion with the 
same potential strength used in the above cited papers. The key point is that the physics 
of the NN central repulsion must be the "same" , with the obvious difference of spin-flavour 
overlaps, as the physics of the K N exotic scattering so that we must be able to derive a 
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common origin for these two experimental facts. We were also able to derive reasonable 
7riV and 7rA couplings for the very same strength of the quark potential which gave good 
KN phase shifts, and all this evidence putted together allow us to hope to obtain, at 
least at a qualitative level, a unified picture for low energy hadronic reactions including a 
microscopically derived NN potential in terms of one single parameter, if the quark masses 
are neglected. 

In section II we briefly describe the model we will be using. The masses of the deltas and 
nucleons, together with the respective Salpeter amplitudes will be obtained in section III. In 
section IV we discuss the set up of the RGM equations appropriated for the NN scattering. 
Section V is devoted to the evaluation of particular overlaps needed to solve the NN RGM 
equation and finally in section VI we present and discuss the results. We conclude with 
some final remarks. 

II. QUARK MODEL WITH CHIRAL SYMMETRY BREAKING 

A. Generator of the Condensed Vacuum 

In this model we assume a new vacuum |0) generated from the trivial vacuum |0), through 
the condensation of qq pairs 

|0) = e Q °~ Qo |0> (1) 

where the generator of the condensed vacuum can be written as 

Qj($) = / d'p^p) M aa ,(9,<f>) b\ cs (p)d\ csl (-p). (2) 

and S are the creation and annihilation operators of q and q, respectively. is the 

radial term of the qq wave function and 

M aiS2 (6,(j>) = -Vto E (lm L lm s \m)()-s l l -s 2 \lm s )Y lmL {6,<P). (3) 

m L ,m a ^ ^ 

stands for the angular term coupled with the spin. The vacuum can be understood in terms 
of a coherent superposition of pairs qq with the vacuum J PC quantum numbers ++ and 
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thus can only have s Lj equal to 3 Po- This coupling is contained in (||). The condensed 
vacuum |0) obtained is orthogonal to |0) if differs from zero. 

The Fock space quark and antiquark annihilators b and d, carry indices for flavour, spin 
and color. In the new vacuum |0), they are related to the old ones by a Valatin-Bogoliubov 
transformation conserving the usual commutation relations, 

{ b \ sc (p), &/vc(p')} = {4*c(p)> df's'c(p')} = 5//'5 ss '<MP - P')- ( 4 ) 



B. Rotation of Spinors 



The field operator ^/ c (x) is defined as 

d 3 p 



^/ c (x) 



(2tt) 3 / 2 



Us(p)b fsc (p) + v s (p)dl sc (-p 



i p.x 



(5) 



Summation over repeated indices is assumed. In order to satisfy the formal invariance of 
ip under the Valatin-Bogoliubov transformation, the rotation of operators in the Fock-space 
must be compensated with an inverse rotation of spinors u and v and we obtain 



MP) 
v s (p) 



1 
1 

V2 



1 + sin (p(p) + J 1 — sin (f(p) p. a 



1 + sin Lp(p) — J 1 — sin (p{p) p. a 



(6) 



where we define the chiral angle 



ni 



(f(p) = arctan — + 2&(p) 



(7) 



With this expression, where m is the current mass of the quark, it is easy to verify that 
for massless quarks the condensation of vacuum generates, through the Valatin-Bogoliubov 



transformation, an effective mass. 



C. Normal Ordered Hamiltonian 



Once the vacuum has been rotated, it is convenient to normal order the second quanti- 
zation Hamiltonian of the model, which is 
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H = J d 3 x[H (x) + # 7 (x)], (8) 

where H is the Hamiltonian density of the Dirac field, 

if„(x) = V f (x) (m,/3 - ia.V) ^(x), (9) 

and an effective instantaneous interaction term 

tf 7 (x) = lJd 3 y ^(x) 7 °y^(x) \/(x - y) ^(y) 7 °y^(y). (10) 

The A a 's are the Gell-Mann color matrices and m q is the quark mass, ip carries indices for 
flavour and color. Flavour is conserved in both H and Hj. 

The normal ordering relatively to the new vacuum gives us the Hamiltonian in the form 

H = H + H 2 + H? + # 4 (11) 

where the subscript is the number of operators in each term. H represents the energy of 
the new vacuum. 

H = N f W J -0^[2C(p)-3(E(p)+m q sin ip(p) + p cos ip(p))\ (12) 
where Nf is the number of flavours, 

C(p) = Jd 3 p'V(p-p'), 

E(p) = sin (p{p) A(p) + cos (p(p) B(p), 
2 



A(p) = m q + -Jd 3 k V(p - k) smip(p), 
B(p) = p + - J d 3 k V(p - k) cos (p(p) p.k (13) 



It is possible to show that this energy turns out to be smaller than the energy of the 
trivial "empty" vacuum (|0); where if — 0). The operator H 2 is given by 

H 2 = J d 3 p E(p) [b\ sc {p) b fsc (p) + 4 sc (-p) d fsc (-p)} (14) 



, where E(p) is the energy of a free quark. 

The anomalous Bogoliubov term has the form 

H2 = J d 3 p [A(p) sin cp(p) — B{p) cos <p(p)) 

■ 'M ss ,b\ sc (i>) 4 s , c (-p) + M* ss ,d fsc (-p) d fs r c (p)] . (15) 

It has the form of the vacuum generator Qq an can creates qq pairs on the vacuum. It is 
now important to remember that the quiral angle, which define the physical vacuum, has 
not yet be chosen. To achieve this we can say that <p(p) minimize the energy of the new 
vacuum. With this condition we get 

A(p) cos ip(p) — B(p) sin<£>(p) = (16) 

and simultaneously we get the anomalous Bogoliubov to vanish. Expression (^) is called 
mass- gap equation and defines the chiral angle ip(p) with an obvious functional dependence 
on the chosen potential. An alternative approach |TT], at least at the BCS level, is to use 
the Ward identity for the H4 kernel in the derivation of the mass-gap equation. 
The term H4 

H, = \J d'pdVq V(q) ) £ : ^ 2 (p, p + q) 0^ 4 (k, k - q) : (17) 

gives 16 terms grouped in ten different four-quark amplitudes. Such amplitudes are combi- 
nations of the following four vertices: 

©c' c (P>P') = "I'(P'K(P) b\ s , c ,(p')b fsc (p), 

e c 2 , c (p,p') = -4(p'K(p) 4 sc (-p)d/ s v(-p'), 

©L(P>P') = «I'(p'H(p) d \s>c'{p')d\sc{-p)> 

e^(p, P ') = ^(p'K(p) d) s , cl {-v')b f sc{v). (is) 
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D. Spontaneous Chiral Symmetry Breaking 



It was shown H that chiral symmetry is spontaneously broken in this model. We briefly 
explain the idea. For massless quarks, the Hamiltonian (|8|) is invariant under the global 
chiral transformation 

if, -> e~ iaaTa ^ ij (19) 

where T a are the SU(Nf) generators and the corresponding axial charges are, 

Ql = J d 3 x^^ 5 T a ^. (20) 

The Bogoliubov transformation on the Fock space of the creation and annihilation op- 
erators can be cast as a rotation through an angle $ = <p/2 and this fact, together with 
the invariance of the field operator ip under these transformations produces a corresponding 
counter-rotation on the spinorial basis. This is enough for the to acquire, for non van- 
ishing <p(p), an anomalous term which is the composite operator for the pion creation and 
thus furnishing the only nonzero contribution for Qg|0). At the BCS level, (Js|0) 7^ implies 
that the pion mass must be zero and the solution of the corresponding Salpeter equation 
for the pion confirms this fact. Furthermore we obtain, as expected, that the pion relative 
wave function is given by sin ip (p). 



E. Harmonic Potential 



In what concerns the phenomenological potential V, both the harmonic and linear confin- 
ing potentials have been studied ||11|| , together with different combinations of Dirac vertices. 



In the present work we will use the Coulomb harmonic confining potential, 

V(x - y) = -\Kl (x - y) 2 , V (q) = ? Kl A q <5(q) (21) 
which furnished the best results in previous calculations. 
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Solving Eqs. (|13|) and fll6|) both with this potential, we can obtain the explicit expression 
for the energy of a free quark, 

_ . , (U)') 2 COS 2 09 

E k = m q sm ip + k cosip — - (22) 

as well as the mass-gap equation, 

p 09 + 2po9 + 2m g p cos 09 — 2p sin 09 + sin2o9 = 0, (23) 

where units of K = 1 are used from now on. 

The mass-gap equation ( |23"|) was solved in Ref [HH, with the boundary conditions 

09(0) = |, O9(p^oo) = 0. (24) 

The chiral angle has been studied in detail but we point out that in the limit of zero current 
quark mass and zero potential, which forces the chiral angle to be zero, expression (|^) yields 
the massless Dirac spinors. The other trivial limit arises for very massive fermions with 
09 = ~, for all small momenta. In this case and for this region of momentum, Eq. (^|) 
becomes an identity. 

III. NUCLEON AND A BOUND-STATES 

Bound-state equation for mesons and baryons were already studied in Ref. p|,|B|,|T2"|j . Here 
we report on more precise results for the hadronic radii and masses, both for the nucleon 
and the A. 

A. The Nucleon Wave Function 

The nucleon internal wave function 0A?(pi, P2, P3), where p^ is the momentum of quark 
i, includes the internal degrees of freedom of the three quarks as color, spin and isospin. In 
the moment space, it is usual to use the ground state harmonic oscillator wave function 

/ 1 x3/4 _ p 2 

^ooo(Pi) = J e ^ (25) 
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as a variational trial wave function for the baryon Salpeter equation, with a variational 
parameter a. It was already shown that it provides an excellent approximation to the 
solution of Eq. (|31]) both for the nucleon and the A. 

Both internal wave functions are orbitally symmetric, SU(A) spin-isospin symmetric and 
SU (3) color singlets. For instance, the wave function of a nucleon is 

4>N = 0ooo(pi) 0ooo(P2) 0ooo(P3) —f=(D s D f + F s Ff) e (26) 

where D (F) is a three quark state which is symmetric (antisymmetric) in the labels of the 
first and second quarks: 

^I = 4|Un + uT-2m) 

F l = -^ (ITT -TIT) (27) 

and Df and Ff are equivalent expressions provided we make the substitutions |— > u, |— > d 
in the case of protons and \— > u, |— > — d in the case of the neutrons. The usual normalized 
color component is totally antisymmetric, 

e = 4g em \ijk). (28) 



B. Baryons Bound State Equation 

For the present study the only quantity needed is the radius of the nucleon in order to 
solve NN scattering. Therefore we solve the bound state equation neglecting the negative- 
energy component, since it corresponds to channels of much higher energy. 

Since our Hamiltonian is instantaneous and there are no negative-energy channels, the 
bound-state equation can be written as 

H\1>) = M\1>), (29) 
where \ip) is an eigenstate of the Hamiltonian with mass M. 
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As in Ref. 0,^,0 we will restrict our calculations to bare baryons i.e. hadrons without 
coupled hadronic channels, in the simplest S-wave configuration, with quarks of the same 
current mass, leading to only one chiral angle (p, the same dispersion relation E(p) for the 
three quarks and only one set of spinors u and v in the vertices 0. 

We begin with the creation operator for baryons which has the form 

*B= / rf 9 p^(Pl+P2 + P3)^(Pl,P2,P3) 

& / lSlCl (Pl) & / 2S2C2 (P2) & L 3C3 (Ps) (3°) 

where ip, \ and e represent, respectively the momentum, the spin-flavour and the color 
components of the baryon wave function. Calculating (H2 + H±)^ B |0) we have for the 
baryon bound state equation: 

[M-3 J B(p 1 )] XsiS2S3 ^(Pi,P2,P3) =-2/ d 3 qV(q) 

[««(pi) w«s(pi -q) M i 2 (P2)^ 4 (P2 + q)]x S 4 S 5s 6 V ; (pi -q,P2 + q,ps) (31) 

with Pi + P2 + P3 = 0, for any flavour indices in the function \. The factor —2 includes a 
contribution of —2/3 from the contraction of the color matrices times the number of quark 
pairs. Thus, the bound-state equation (|3"T|) has the form of a Schrodinger equation with an 



effective interaction 
2 r 

V = -- J d 3 qV(q) [?4(pi) M S3 (p! - q) wj 2 (p 2 ) w S4 (p 2 + q)] Xs 4 s 5 s 6 ^(Pi ~ q, P2 + q,Ps)- 

(32) 

Calculations made with the harmonic potential ( pT| ) give 

= v!L + v*i + + v£g (33) 

including the spin independent, hyperfine, tensor and spin-orbit interactions [O . 
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C. Nucleon and A masses 



In order to check the consistence of our model, we calculate the N and A masses. In 
the approximation of N and A with only one S-wave channel, the spin-orbit and tensorial 
terms vanish and the remainder bound-state equation is 



{3^,)-M-^„ + g (y w + 3|1 ' s ;r (f " )l 



+ 



l-\s(s + i) 



[1 - sin (pfa)] [1 - sin if(p 2 )] \ ^(pi, p 2 , p 3 ) = 0, (34) 

P1P2 J 



where S is the total spin and Vp 12 stands for the Laplacian with respect to the relative 
momentum (pi — pa)/2. The mass of the baryon is then given by 

. , « , , 3 ~ 3. ,. 9 . . 3 fl — sin <p(pi)l . 

M = (3E( Pl ) — M — -Vj 12 + -(^) 2 (Pi) + -2 ;J ) Pl+P2+ P3=o 

z 4 p 1 



l~l S{S + 1] 



([1 - sin<^(pi)] [1 - sm(f{p 2 )} ^^)p 1 +p 2 +p 3 =o (35) 

P1P2 



where the following notation has been introduced: 

(/(Pl,P2,P3))p 1 +p 2 +p 3 =0 = J d 9 ?)5{Vi + p 2 + p 3 ) /(Pl,P2,Ps) l^(Pl,P2,P3)| 2 (36) 

Since we do not have an analytic expression for the chiral angle <p, the calculations were 
done with the good approximation 

sin^(p) = e - p2 ^ (37) 

where p = 0.582. 



Figure |VII| shows the results for the NA masses, in units of Kq, as a function of the 
variational parameter a, for the case of massless current quarks. The results are presented 
in table |V1J| considering the harmonic potential with a strength constant 247 MeV. 

The bare nucleon mass obtained gives a mean radius of 0.66 fm, which is reduced when 
the strength of potential is increased. As an example, Kq = 290 MeV gives 0.57 fm. The 
introduction of finite current quark masses almost does not change the results we have 
obtained. 
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The difference between the bare nucleon mass and its physical counterpart is a conse- 
quence of the coupling to Nir channels which is expected to lower the nucleon bare mass 
somewhere between (—300) and (—400) MeV ||13|| . The effect of coupled channels should 
also contributes to the physical N — A difference (100 — 200 MeV). 



IV. RGM EQUATION FOR NN SCATTERING 

It is well known from potential scattering theory that a free spherical wave converg- 
ing towards the origin is almost unaffected by a potentials which range is of the order of 



L(L + l)/k. The repulsive short range potential between two nucleons has almost no ef- 



fect on waves for which J L(L + l)/k > 0.8. This fact motivates the choice of the relative 
S wave to test the repulsive core of NN scattering also revealed by experimental low-energy 
scattering phase shifts in 1 Sq and 3 Si partial waves. 

To study the central microscopic NN repulsion it suffices to consider the hyperfine quark- 
quark potential V^ 2 , 

q ^ 2 2 p iPj K ' 

As we are neglecting the pionic cloud responsible for the long and intermediate range forces, 
it is obvious that we will not be able to reproduce the NN phase shifts. However, the results 
will allow us to show how much the quark structure of the nucleons influences the NN S-wave 
scattering. 

To investigate the Nucleon-Nucleon scattering we make use of the Resonating Group 



Method (RGM) introduced by Wheeler |14j], and profusely applied by several groups 



PJ3| JT5|JT6|1 . In this method the two clusters relative motion wave function is obtained by 
requiring it to minimize energy of the system. Then, we must to minimize the matrix element 

I = (N t N 2 \(H - E)]^ N 2 ) (39) 

where H = H 2 + H± is the Hamiltonian in the second quantization formalism ones, E rep- 
resents the total energy of the system and \N\N 2 ) describes the six quarks system clustered 
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in nucleons Ni and N 2 . 

The Pauli principle is taken account by demanding the antisymmetry of the total six 
quarks ensemble which arises naturally through the Wick ordering of quark fields. Since the 
wave functions of nucleons are already antisymmetrized, we obtain, 

(NiN 2 \(H - S)|JViJV 2 ) = (N 1 N 2 \(H - E) 36(1 - 9P M )(1 - V)\N X N 2 ) 

= (N 1 N 2 \(H — E) A\N 1 N 2 ) (40) 

where Pij stands for the permutation of the quarks i G N\ and j G N 2 and V = -P14-P25-P36 
gives the permutation of the nucleons 1 and 2. A is the usual designation of the antisym- 
metrizer of six quarks in two nucleons. The operator H is defined in first quantization. 

At this stage it is interesting to see that the factor (1 — V) implies, in a natural fashion, 
the usual antisymmetry condition at the baryon level. In fact, (1 — V) acts on the system 
of two nucleons and makes the matrix element vanish unless the condition L + S + T = odd 
holds. 

The Jacobi coordinates adapted to this problem are 

JH-P2), Pa = ^(Pi + P2 - 2p 3 ) 
>4-P5), P 7 = ^t(P4 + P5 -2p 6 ) 
>i + P2 + Ps) - ^(P4 + P5 + Pe) 

Pcm = 7/g(Pl + P2 + P3 + P4 + P5 + Pe)- ( 41 ) 

As a result of this choice, each nucleon is described by one of the two sets {(p, A)}, {(£,7)} 
of coordinates. The relative motion is described through p r . 

The Hamiltonian in first quantization, in the center of mass frame, becomes 

H = ^ + ^ + ^ + ^ + ^ + ^2^3 (42) 

i<j 





1 


Pp = 


72 




1 


Pe = 






1 


Pr = 


7q 



The subscripts p, A, £,7, r and cm correspond to the Jacobi coordinates of equation ( flTD 
and Vij is the potential of interaction between the quarks i and j. Using a more explicit 
notation, 
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/ = 36(x(p r )l(0 1 (6)0 2 (6)l(^ - E){1 - 9P 14 )(1 - ^)|0 1 (6)0 2 (6))lx(p r ))p cm =o (43) 

where the arguments £1 and £ 2 denote the sets of internal degrees of freedom in the wave 
function 1 ' 2 of each nucleon, including P p ,Pa,P$,P7 in the momenta space, and x(Pr) is 
the wave function for the relative motion of the two nucleons. 
Expanding x(Pr) in spherical harmonics, 

X(Pr) = E E Xim(Pr), (44) 
m I 

where Xim{Pr) includes both radial and orbital components, we obtain that for the particular 
case of S relative wave function Xs(Pr) = Xs(— Pr), and we can drop the permutation V 
between the two nucleons, providing the product 0*(£i) 2 (^2) of the two wave functions 
is antisymmetrized. The RGM equation is obtained minimizing the matrix element / with 
respect to Xi 
8T 

^ = o^ (0 1 (6)0 2 (6)I(^-^)(1-9P 14 )|[0 1 (6)0 2 (6)])|X5) = o. (45) 

It is easy to verify that the specific dependence of color in yields a null matrix element 
when the permutation is not present. The nucleon mass Mjv is defined by 

E N \N X ) = (T p + T A + V 12 + V 13 + V 23 )\N) = M N \N 2 ), (46) 

with an analogous expression for the second nucleon. 

Now we introduce the eigenfunctions 0" Zm = <p a of the harmonic oscillator with parameter 
a as a complete set for the expansion of the relative wave function Ixs) 

lxs> = EI&W«H&Wto>. ( 4? ) 

ab 

Using the fact that E r = E — 2 M N , we get 

[(Tr-ErXl + a^&bfi+VabltaHM \Xs) = 0, (48) 

where 
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a ab = -9(0 a ( Pr )|(0 1 (6)0 2 (6)|i 3 i4|[0 1 (6),0 2 (6)])l06(p r )) = -9(P 14 ), 

V ab = -HMPrWi&A&l E ^ i P 14 |[0 1 (ei),0 2 (6)])|06(Pr)>. (49) 

i€Ni,jeN 2 

We have to note that Eq. (fi?D, associated to expansion (pB[), implies that = n ,o,o- We 
will retain a = b = (0, 0, 0) as an approximation that is reasonable due to the fact that the 
overlap kernel being compact, happens to filter out higher wave excitations quite effectively. 
Now, if we use 

|^) = (1 + ( too|0ooo)(0ooo|)|x) (50) 

where c^ooo is the wave function of the tridimensional harmonic oscillator with the relative 
Jacobi momenta p r it is possible to rewrite Eq. (f|^) as a non-local Schrodinger type equation: 

f Tr - Er + -^^|0 OOO )(0 OOO |) |V>) = (51) 
V 1 + o"oo / 

Of course, wave function fl5"0|) gives the same phase shifts than \x)- 



V. EVALUATION OF a AND V 

The nine diagrams we have to consider are given in Fig. |V11| . The values of <Too for 
the color and spin-flavour spaces were directly calculated numerically, performing all the 
permutations involved. To obtain the value of t>oo i n the same spaces, we must consider the 



nine graphs of Fig. [VI I| . In order to simplify calculations, the interactions of color and spin 



can be commuted in permutations given by relations: 



A? A? 1_„„,^ _„ 1, 



' 'i 3 _ _( pcolor I or _ 2. 
2 • 2 - 2 ^ij +*>ij y 

U 3 = \{Pr-\) (52) 
The evaluations in momenta space could be made with the appropriate graphical rules 



17| . But in our order of approximation, we do not need to use them. For instance, to obtain 



(Pu), we just have to note that P14 does not change the momenta part of the wave function 
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which describes six quarks grouped in two nucleons and with relative wave function 0ooo- 
This result is instantaneously obtained by inspection of expression 

{4>Nl4>N2)mom = 0OOo(Pp) 0OOo(Pa)0OOo(P?) 0000 (P 7 ) 0000 (Pr) (53) 

Therefore a — 1. 

For each one of the nine terms of vqq, the procedure is the same we describe. We start 
observing that for the reason given above, permutation P14 can be removed. In momenta 
space, each term of t>oo reveals to be proportional to 

Voo = (0ooo(p,)|(0 1 (6)0 2 (6)l(^) m o m |0 1 (ei)0 2 (6))l0ooo(Pr))p cm =o, (54) 

where (V^) mom is the momentum dependent part of the hyperfme quark potential. 

A numerical evaluation gives Voo — — 0.0671Aq- The calculations were done for the cases 
(T, S) = (1, 0) and (T, S) = (0, 1) and the results are in table [VI]] , where the evaluations in 
the spin-flavour space must be multiplied by a factor (1/18) 2 . 

For instance, we have, for (T,S) = (1,0): 

^00 — (Pu)c-(Pl4)sXf-(Pu)mom. = r 

y 

5 

V00 = -9 ((Vij Pli) c- (Vij P U ) sx f-(VijPu) mom.) = --VOO (55) 



These results are summarized in table VII. 



VI. RESULTS AND DISCUSSION 

The value of K , the unique parameter of this model, was fixed to K = 247 MeV in 
order to obtain good KN phase shifts. This value also gave good f^N and /tvttA and was 
used successfully by the Orsay quark group . We obtain a = 296.4 MeV and —Voo — 16.6 
MeV. In the most conventional quark models, the oscillator parameter has values of 250 
MeV JL]| to 410 - 420 MeV 1I6|J2(J. 

Since the potential is separable, the phase shifts were obtained analytically 
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5(z) = -atan- 2 ^^ex P (-z>) 

with 

, = * E r =M (57) 
1 + o"oo Oi 2m 

where m is the physical reduced mass of the system of two nucleons and E r stands for 

the energy of the relative motion in the center of mass frame. The factor V6/2 is needed 

because p r is a Jacobi coordinate. The phase shifts obtained are presented in Figs. |V11| for 

the (T, S) = (1, 0); (0, 1) channels. 

In order to derive an equivalent local potential Uelp defined, according to Ref. ||, as 

2/i 

the nonlocal Schrodinger equation (|5lD was solved. This local potential happens to be 
energy-dependent. The numerical accuracy was tested by reproducing numerically the pre- 
vious phase shifts. 

It is expected in a microscopic calculation of the AW interaction, to obtain for the short 
range part a non-local description, besides the usual meson-exchange picture which is also 
natural in this framework. There are a few realistic potentials that fit in this description, 
with a different height for the core, describing very well the two body data. The latter does 
not seem to be very sensitive to the height of the core. The Nijmegen potential was chosen 
for comparison because it has a version, the Nijml |Tj|, with a very soft core, comparable to 



^_ V 2 + E r - U ELP (r) 



V>(r) = , (58) 



ours, and has an optimal x 2 P er datum for np data. In Fig [VI I| , both our equivalent local 
potential Uelp and the Nijml potential are depicted. It is well known that they should differ 
in the medium range, due to the meson exchange contributions which are present in any 
realistic potential, and that we have not yet included. Nevertheless it is remarkable that, 
contrary to past experiences, such a soft repulsive core can be accomodated to reproduce 
experimental phase shifts. 

Almost all the AW potentials obtained in the framework of quark models have compara- 
ble range values, but they differ largely in the magnitudes. For example, Faessler et al and 
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Oka P,|T6l| find potentials with peak values between 600 and 1300 MeV. Ribeiro @ built a 
square potential of magnitude 500 MeV that was approximately phase shift equivalent to 
the ones obtained solving the RGM equation. The Liberman potential has a peak of 450 
MeV. Of course, a direct comparison of the potentials is not absolutely correct because the 
methods used to derive each potential are quite differents. In order to avoid this problem, 
we can compare the equivalent local potentials obtained by Suzuki and Hecht [ffiq] , who used 
a common method for the models we have referred above. However we can accept this large 
range of results if we take in account the small relative weight of the potential near the 
origin. 

Although it is usually expected that quark exchange processes give the most relevant 
contribution to the NN repulsion, u and p can still give contributions to the short range 
potential. These mesons, and the other possible ones required to describe the intermediate 
and long range attraction of the nuclear force, can be introduced in a consistent way within 
this model. 



VII. CONCLUSION 

In this work we have derived, from a chiral invariant quark model, the short range part 
of the NN interaction. There is only one free parameter in this model, which was already 
fixed in previous calculations. A non-local potential was then derived, with a soft core, 
consistent with the recent Nijml realistic interaction, that fully describe the experimental 
NN data. This is part of a wider program to describe the full nucleon-nucleon interaction 
within a quark model. 

Finally we find remarkable that we are able to understand, with the same parameter, 
the NN repulsion and, at the same time, to reproduce the meson and baryon spectra, the 
vector meson decays, the f n-kN coupling and to obtain a reasonable KN repulsion. 
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FIGURES 

FIG. 1. Masses of N and A as a function of the variational parameter a. Both are in units of 
K = l. 

FIG. 2. The nine diagrams that contribute to the potential overlap. 

FIG. 3. Phase shifts of NN scattering for (T, S) = (1, 0) and (T, S) = (0, 1) channels. 

FIG. 4. Our equivalent local potential (solid line) and the realistic Nijml potential (dotted 
line). The deviation for large distances is due, as usual, to attractive meson exchange. 
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TABLES 





N 


A 


Mass (K ) 


4.88 (1207 MeV) 


5.08 (1255 MeV) 


a (K ) 


1.20 


1.21 



TABLE I. A, A masses and inverse radius a in units of Kq = 1. Values in MeV were obtained 
with A =247 MeV. 





(T,S)=(1,0) 


(T,S)=(0,1) 




color 


spin-flavour 


mom. 


color 


spin-flavour 


mom. 


(Pu) 


1 
3 


-12 


1 


1 
3 


-12 


1 


{VuPu) 


4 
9 


93 


Voo 


4 
9 


57 


Voo 


(Vi 5 Pu) 


2 
9 


-21 


Voo 


2 
9 


-21 


Voo 




2 
9 


-21 


Voo 


2 
9 


-21 


Voo 


(V24P14) 


2 
9 


-21 


Voo 


2 
9 


-21 


Voo 


(^25^14) 


1 

9 





Voo 


1 

9 


5 


Voo 




1 

9 





Voo 


1 

9 


5 


Voo 


(L34P14) 


2 
9 


-21 


Voo 


2 
9 


-21 


Voo 


(V 35 Pu) 


1 

9 





Voo 


1 

9 


5 


Voo 




1 

9 





Voo 


1 

9 


5 


Voo 



TABLE II. Contribution of each diagram in color, spin-flavour and momenta spaces. 
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(T,S) = (1,0) 


(T, 5) = (0, 1) 


O"00 


1 

9 


l 

9 


^00 


-|V o 


104 v 


\ — voo 
A - 1+ctoo 


-fVoo 


-§V 00 



TABLE III. Final values for <7oo and vqq. The relevant quantity to obtain phase shifts is A. 
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